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THE KLEIN-GORDON EQUATION ON THE TORIC 
ADS-SCHWARZSCHILD BLACK HOLE 

JAKE DUNN AND CLAUDE WARNICK 


Abstract. We consider the Klein-Gordon equation on the exterior of the toric anti 
de-Sitter Schwarzschild black hole with Dirichlet, Neumann and Robin boundary con¬ 
ditions at X. We define a non-degenerate energy for the equation which controls the 
renormalised norm of the field. We then establish both decay and integrated de¬ 
cay of this energy through vector field methods. Finally we demonstrate the necessity 
of ‘losing a derivative’ in the integrated energy estimate through the construction of 
a Gaussian beam staying in the exterior of the event horizon for arbitrary long co¬ 
ordinate time. 


1. Introduction 

Among asymptotically flat spacetimes satisfying the vacuum Einstein equations: 

( 1 ) Ricg = 0 , 

the black hole solutions occupy a privileged position. It is conjectured that the family 
of rotating black holes described by the Kerr metric represent the final state of gravita¬ 
tional collapse. Within the class of spherically symmetric solutions with scalar matter 
it is possible to establish this fact rigorously [T]. Establishing such a result without an 
assumption of symmetry appears a very challenging task. Less ambitiously, one might 
hope to show that the Kerr family of black holes are in some suitable sense stable as 
solutions of Q against small perturbations to the initial data. In the absence of a black 
hole, the nonlinear stability of the Minkowski spacetime against small perturbations was 
established in the work of Christodoulou-Klainerman [2]. For black hole spacetimes, 
the nonlinear stability against generic small perturbations remains an important open 
problem, although see lam for an alternative approach. 

In the linear setting, considerable progress towards establishing the stability of the 
black hole spacetimes has been made over the last few years, culminating in the recent 
proof of decay for solutions of the scalar wave equation on any fixed subextremal Kerr 
black hole background [5]. Even this highly simplified problem nevertheless requires 
a very careful analysis. The key issues to be understood are those of trapping and 
superradiance. The issue of trapping is related to the presence of null geodesics which 
orbit the central black hole, neither escaping to infinity, nor falling through the black 
hole horizon. Trapped null geodesics are an obstacle to linear decay, resulting in decay 
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estimates that ‘lose derivatives’. The other problem to be confronted is the superradiant 
effect, which manifests itself in the absence of a globally timelike Killing field. As 
a consequence, the natural conserved ‘energy’ is not coercive and does not control the 
solution. Overcoming these issues requires a rather subtle argument involving the nature 
of the set of trapped geodesics and its interplay with the superradiance phenomenon. 

In this paper, we shall initiate the mathematical study of the Klein-Gordon equation 
on a family of black hole backgrounds which, largely, evade the complications of trapping 
and superradiance present for the Kerr family of black holes: the planar, or toroidal, 
AdS-Schwarzschild black holes. These are solutions of the vacuum Einstein equations 
with a negative cosmological constant: 

(2) Ricg = hg, A < 0. 

The presence of a negative cosmological constant for our purposes has two main effects. 
Firstly, the character of null infinity changes: it becomes a timelike surface, meaning that 
it is necessary to impose boundary conditions in order for the Klein-Gordon equation to 
define a well posed evolution problem. Secondly, for A < 0 it is no longer the case that 
the horizon of a compact black hole in four dimensions must be spherical: versions of the 
AdS-Schwarzschild black holes exist whose horizon is a surface of arbitrary genus. We 
shall consider the case of a toroidal AdS-Schwarzschild black hole, with isometry group 
M X T^. The spherical case has been studied in PE!. 

The main difference between the toric AdS-Schwarzshild black holes and the spheri¬ 
cal AdS-Schwarzschild black hole can be seen at the level of the null geodesics. In the 
equation governing null geodesics in the toric black hole background, there is no ana¬ 
logue of the centrifugal force that appears for the spherical black hole. Grudely, there 
exist trapped null geodesics in the spherically symmetric case owing their existence to a 
combination of the centrifugal repulsion and the confining nature of the AdS boundary. 
For the toric black hole, the centrifugal force is absent, and so no trapped null geodesics 
exist. This has a profound effect on the behaviour of solutions to the Klein-Gordon 
equation on this background. 

Aside from the fact that solutions of (§ are of interest in the context of classical 
relativity, they have also received considerable attention in the context of the putative 
AdS/GFT correspondence. This conjectures a relation between gravitational theories 
with negative cosmological constant, and conformal field theories in one fewer dimension. 
In this context linear fields on the planai]^ AdS-Schwarzschild solution are studied, as 
they are conjectured to be dual to conformal field theories in Minkowski space at finite 
temperature. 

The main results. We shall first consider solutions of the Klein-Gordon equation: 

. . oi 

(3) DgU-I-= 0, 

where a < |. (When dealing with integrated decay due to a technical limitation we will 
restrict to a* < a < | where a* is a constant with the property a* G (|, |g), {a* 1.46). 

Our first result is an extension of the results of [8| to the toric AdS-Schwarzschild black 

^The difference between planar and toric AdS/Schwarzschild is purely topological: the toric solution 
arises by periodically identifying the planes of symmetry present in the planar solution. 
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hole. That is to say, for a stationary time foliation which is regular at the horizon, we 
shall establish the following result: 


Theorem 4.1 


Suppose u is a solution to (§ satisfying suitable (Dirichlet, Neumann 
or Robin) boundary conditions at infinity. Let £[u] be the renormalised energy density 
of the field u. Then for any T > 0 we have: 


l{t=T} 


£[ 


U\ 


< 



e[u], 


with the implieit constant independent ofT. 


Here and elsewhere, £\u\ is an energy density which does not degenerate at the event 
horizon, so we gain full control of all derivatives up to, and including, the horizon. 
Having established boundedness, we then consider the question of decay for solutions of 
Q. We establish an integrated decay estimate 


Theorem 5.1 . Suppose u is a smooth solution to (i satisfying Diriehlet, Neumann 


or Robin boundary conditions at infinity. Let £[u] be the renormalised energy density of 
the field u. Then for any T > 0 we have: 


f ^ 

J{o<t<T} 



£[u\, 


with the implieit constant independent ofT. 


Here r is a radial coordinate with null infinity corresponding to r —oo. Clearly, the 
weight on the left hand side of this estimate is weaker than that on the right. In fact, 
this degeneration occurs only for derivatives tangent to X. We can improve the weight 
at the cost of losing a derivative: 


Theorem 5.2 . Suppose u is a smooth solution to ([^ satisfying the same boundary 
eonditions as above. Let £[u\ be the renormalised energy density of the field u. Then for 
any T > 0 we have: 


£[u] 


< 


J{0<t<T} J{t=0} 

with the implieit constant independent ofT. 


{£[u] +£[ut ]), 


It is likely that one does not need to lose a whole derivative to improve the weight, 
but we shall not pursue this point here. 

In order to establish the results above, we combine the vector field method with the 
renormalisation methods of [9]. The boundedness proof follows a very similar approach 
to that of |8]. To establish decay, we make use of currents constructed from Morawetz 
vector fields. In contrast to the spherical black holes, the absence of a photon sphere 
simplifies the construction of appropriate currents. 

Combining the integrated decay estimate of Theorem 5.2 with the result of Theorem 
and an argument based on the redshift, we can establish: 


4.1 
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Theorem 5.3 . Suppose u is a solution to ([^ satisfying suitable (Dirichlet, Neumann 


or Robin) boundary conditions at infinity. Let £[u] be the renormalised energy density 
of the field u. Then for any T > 0 we have: 


l{t=T} 


£[u] 


< 


1 


(l + T)' 


E 

k=0 


/ ^ 


J{t=0} 



with the implieit constant independent ofT. 


It is worth contrasting this result with that of PE] for the spherical AdS-Schwarzschild 
and Kerr AdS black holes. They established that in this background solutions of the 
Klein-Gordon equation satisfy a decay estimate of the form: 

and moreover this cannot be improved. The slower rate of decay here is due to the exis¬ 
tence of null geodesics which never cross the horizon: they are stably trapped between 
the centrifugal barrier and null infinity (which for these purposes may be thought of as 
a reflecting barrier). 

The fact that our estimates exhibit a loss, either in weight or in derivatives, is due to 
a ‘trapping at infinity’. There exist null geodesics which start far from the black hole, 
moving orthogonally to the radial direction, and which take an arbitrarily long time to 
cross the horizon. Using this fact, one can establish using the Gaussian beam methods 
of [inj that: 

Theorem |5.4| . There exists no constant C > 0, independent ofT, such that the estimate 

[ £[u\ ^ C f £{v\, 

J{0<t<T} J{t=0} 

holds for all smooth solutions u of (§. 

This result is very robust, in particular it does not depend on the boundary conditions, 
nor on the value of the Klein-Gordon mass. 


Outline of the paper. This paper consists of four main sections. The first one defines 
the spacetime, sets up the appropriate hypersurfaces we will need and states the relevant 
divergence theorem. The second section sets up the initial boundary value problems 
(IBVP) we are interested in, the appropriate re-normalisation scheme for well posedness 
as seen in [9| and defines the necessary tensorial quantities for proving decay. The third 
section proves energy decay for the IB VPs, taking the approach of [8|. The final section 
then seeks to obtain qualitative time decay rates for the energy of solutions of the IB VP 
in the time coordinate using vector field methods found in m and [T^ . 


2. Toric AdS Schwarzschild black hole 

2.1. The Manifold. We define the exterior of the toric AdS-Schwarzchild black hole 
with mass M and AdS radius I to be the following manifold with boundary 

Nt = Mj>g X X 
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with Lorentzian metric 


-2M r2 


g = - 


+ -j2 ) dt'^ + 


4Mt^ 


■dtdr + 


f2Ml^ P\ ,2 2 , 2 ^ 

+dy )• 


Where denotes the two dimensional torus and r+ = {2MP)3 is the event horizon for 
this space. This metric was first stated in the papers [ElIIl!. 

In these coordinates {t, r, x, y) we make the usual periodic identification for the torus 


X ~ X + pi, y ~ y + P2- 

As the toric domain is compact integrating over it shall prove no problem and these 
periods will typically not appear explicitly. There appears to be no obstacle in extending 
our results to the planar black hole (without identification of x, y), provided one assumes 
some decay for the field in the x, y directions. 

For later asymptotic analysis it is also helpful to know the expression of the cometric 



2.2. Hypersurfaces and Measures. In this paper we will make extensive use of the 
divergence theorem. With this in mind it is helpful to establish a few hypersurfaces and 
measures. 

The volume element for this space is 

d'^V ol = r^drdtdxdy = r^drj. 


We introduce the following spacetime slab 

~ '^Ti<t<T2 X I^r>r+ X T^. 

We now define the hypersurfaces we will need 

• Si the hypersurface of constant t. This surface has future directed unit normal 
given by 


n 


= - 


<^rt 




--dr 


2l^M P ^ 

-TT- + -^dt - 


21M 


y/2PMr + 


dr 


and induced surface measure 

dSs^ = -s/—g^^r'^drdxdy 


= iJ + r^'^drdxdy. 


A simple calculation shows that Si is a regular spacelike hypersurface up to and 
including the horizon. The following notation will also be useful 

[.Rl A2] 


Si n {Ai <r< R 2 }. 
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• Tij. the hypersurface of constant r. This surface has unit normal given by 



and induced surface measure 

dSs^ = \/W^r‘^dtdxdy 

drdxdy. 

Notice that m becomes singular and degenerates as we approach the hori¬ 
zon. The combination m^dS-£^ is well behaved and gives the appropriate normal 
volume element to the surface. We again will make use of the notation 

^[Ti,T2] = n {Ti < t < Ts}. 

• The surface denotes the intersection of St and S^. It has induced surface 
measure 

dSrp 2 ^ = r^dxdy. 

• We will also denote the hypersurface of the event horizon by 

n = {r = r+}, 

and define null infinity X formally as 

X = {r = oo}. 

While working in this spacetime it is helpful to initially restrict to a finite region of the t 
and r coordinates then take limits to recover the full space. With this in mind we define 

B = {{t,r,x,y) G [Ti,T 2 ] x [Ri,R 2 ] x T^}. 


I 


Figure 1. Penrose diagram of the spacetime 
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We now state the divergence theorem for this setting. If we take to be a vector 
field on M then 


(4) 


'B 




>-^ 2 ] 

"^2 




n[«l.«2l 

-Ti 




+ f Jn'm'^dS'Srf — [ Jnm^dSY,r, 

fll R 2 


Providing the limits exist we can extend this to the spacetime slab ■M.\yX\,T 2 \ by sending 
Ri —> r+ and R 2 —)• 00 . 


To simplify our dealings with tangential terms we will denote the connection of the 
induced metric on tori of constant t and r by y. We note that in our coordinate system 
for a function u G C^(AI) we have 

^ ['^l + ^D ■ 


3. The Klein-Gordon Equation 

3.1. Klein-Gordon Eqnation. The Klein-Gordon for an asymptotically Anti de-Sitter 
spacetime is given as 

. . o 

(5) □gU-h^?X = 0, 

where a < | obeys the Breitenlohner-Freedman bound |15] . 

Using the coordinate expansion 


UgU = 


d, 


VW\ " 

for the wave operator, Q takes the form: 

/ 2ml‘^ P\ 1 / o / —2m 

— I -=-1- n I Utt+^Or 




u 




r2 


u. 




2mf 


a 


Uft 


^t + :p2^{x,y)U+j^U = 0 


We note that second order radial derivatives degenerate as r —)> r_|_. If we express this 
PDE in the form 


-utt + But + Art = 0, 


for some spatial differential operators L of degree one and two respectively, we find 
L is not strongly elliptic at the horizon and that the standard energy methods are 
insufficient to prove boundedness of the full norm. We overcome this by exploiting 
the redshift effect for black holes. The details on how this can be done for general black 
holes can be found in m- We must also confront the issue that standard energy fluxes 
of this PDE diverge as r —>• 00 this is fixed by a renormalization process, in particular 
by reformulating the problem in terms of the twisted derivative. 







3.2. The Twisted derivative. Due to the fact that asymptotically anti de-Sitter spaces 
do not admit a Cauchy hypersurface, in order to have a well-posed problem it is necessary 
to prescribe boundary data on X. While the standard dt energy currents are sufficient to 
establish well-posedness for Dirichlet type data, for more general boundary conditions the 
solution has less radial decay and the standard energy diverges. This was resolved in [9] 
for the range a G (|, |) by treating the well-posedness of the problem in asymptotically 
AdS spaces with the use of a re-normalisation scheme. The core idea is to reformulate 
the energies in terms of ‘twisted’ derivatives 

= /V^ {f-\) , 

for a ‘twisting’ function / > 0 which captures the decay of the field near I. 

Defining as the formal adjoint of V with respect to the spacetime inner product 

we note we can rewrite the Klein-Gordon equation in the form 

-Vu = 0, 

where 

3.3. Boundary conditions. It will be helpful to define K>0bya = | — We 

say that u G obeys Dirichlet, Neumann or Robin boundary conditions if the 

following hold 

• Dirichlet, k > 0 and 

3 

r'2~'^u —)> 0, as r —)> oo, 

• Neumann, k G (0,1) and 

> 0, as r —)• oo, 

• Robin, K G (0,1) and 

r+ (3r'2~'^u —)• 0, as r —oo, 

where /3 G 

We will throughout this paper assume as part of our definition of Robin boundary 
conditions that dtf3 = 0 and /3 > 0. 

3.4. Well-posedness and asymptotics. In this section we state the well-posedness 
results and the asymptotic behaviour of the solution as found in [9j. 

Firstly let S be a smooth spacelike hypersurface which extends to I and meets it or- 
thogonalljj^ We let ns be the future directed unit normal of S and define 

hs = rns, 

then let denote the future Cauchy development of S together with the portion 

of I lying in the future of S. 

O 

with respect to any conformal regularisation of the boundary. 




9 


We choose our twisting function / such that /ra ^ = 1 + 0 (r ^)asr—7-ooin order to 
define the norms 

and the space Hq{Ti, k) as the completion of the smooth functions supported away from 
X. Different choices of twisting function / satisfying the same asymptotic condition give 
rise to equivalent norms. 

Theorem 3.1 (Well-posedness and asymptotics). 

• Let uq G Ml G X^(S). Then there exists a unique u such that m|s = 

uq, usm|s = ui which weakly solves 

in X+(S) with Dirichlet boundary conditions on X. If S is any space like surface 
in X’+(S) meeting X orthogonally then G hj:u\s G Lf{S). 

• Let Uq G Ml G Lf{'E) and 0 < k < 1. Then there exists a unique u 

such that m|s = uq, hsM|s = ui which weakly solves 

+ (i ^ = 0’ 

in X>+(S) with Neumann or Robin boundary conditions (for given (3) onX. If S 
is any space like surface in X’^'(S) meeting X orthogonally then m |5 G II}{S,k), 
hY.u\s G X^(5). 

If the initial conditions satisfy stronger regularity and asymptotic conditions, then 

m |5 G (S), hsu\s G for any integer k > 2 and we obtain an asymptotic 

expansion 

u = X- {«- + 0(r->-)) + -N („+ + 0(r-‘)) . 

Where the functions u~ G H^~^{X), u'^ G H^~‘^{X) satisfy 

u~ = 0 if u satisfies Dirichlet data, 

= 0 if u satisfies Neumann data, 

2ku~^ — [3u~ = 0 if u satisfies Robin data. 

Remark 3.1. We will for the remainder of this paper assume our solutions are smooth, 
and admit asymptotic expansions to all orders. (Such solutions can be constructed from 
sufficiently smooth initial data.) This assumption can later be removed by a density 
argument. We state the asymptotics for the solution of the Dirichlet, Neumann and 
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Robin problems, when k G (0,1); 
Dirichlet: 


u = O [r 2 

5 


VtU = 


\fu\ = O (r- 2 -^) , VrU = 0 (r -2 


u = O yr 2 ' 

\^u\ = O (r~i~ 


Vtu = 0{r-l+^ 


VrU = O [r 2 


Neumann/Robin: 


( 6 ) 


3.5. The Twisted Energy Momentum Tensor. When deriving energy estimates 
for wave equations on the exterior of black holes one typically considers the energy 
momentum tensor. By choosing suitable vector fields as multipliers and commutators 
one can prove boundedness and decay of solutions (for a discussion of the history of this 
approach in the black hole setting see m)- In the case of AdS spaces with Neumann 
and Robin boundary data one finds that for the standard energy estimate obtained by 
using the vector held T = dt as a, multiplier, the energy is no longer hnite. This is 
resolved by the introduction of the twisted energy momentum tensor dehned as 


+ Vu" 


where again; 


V = - 


V/.V^/ 

/ 


a 




Unfortunately the divergence of the twisted energy-momentum tensor (in contrast to 
the untwisted version) is no longer vanishing for solutions to However, it does have 
some useful properties 

Lemma 3.1 (taken from [8j). 

• For/G C^{M) 

= (-vt - Vi!:) 

where 

5 rn Vt(/U) ,2 , Vt/^ , 

Sum = — 2 ^—P + i’- 

• For u a solution to the PDF and X a smooth vector field. Defining 

+ X’^S^u], 

where ^TTf_iu is the deformation tensor 


T^liU — 2 (^/lAjy + X yX^) — ^ (Cxg) 




x^jmu] = k^[u]. 


we have 
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• If f is chosen such that V >0 then satisfies the dominant energy condition, 
i.e. for a future direct causal vector field Y we have that —is future 
directed and timelike. IfY is timelike then T^i^Y^Y'^ controls all the first order 
derivatives ofu. 

As we have that and thus K^[u] depend only on the 1—jet of u we get that J^[u] 
is a compatible current in the sense of Christodoulou HZ!. 

Importantly if we have that Z is a Killing field that preserves /, i.e. Cz{f) = 0 then 
J^[u] is a conserved current. 

We remark that the renormalization by the twisted derivative encompasses the holo¬ 
graphic renormalization of [18] and m- 


4. Energy Decay 

With the notions of the previous sections in mind we are now in a position to prove 
an energy decay result for the Klein-Gordon equation in this setting. This follows the 
same method as in [8] which deals with the spherical case. 

Firstly we choose our twisting function. Experimentation suggests the following function 


/(r) = r 2 +'^. 


This choice gives 


, (3-2k)2m ^ 

“ 2r3 ^ 

and we can easily verify that for the timelike vector field 

T = dt, 

we have 

Cxif) = 0, 

so that 

Jj = 0. 

We now integrate over B to get the following identity 

f Jln^^dS^^ - [ = f - [ j'fm>^dSj:„ , 

/ [Hl.flal ^^2 / [fli.fla] ^^1 LlTi,T2] ^«2 / [Tl.Ta] 

t/ Z-Jrj-i t/ tJ 2—i Ty ^ Z—t D 


«2 


Through a long but straight forward calculation one arrives at 


£t{n- [Ri,R2]): = 




—g^^{'Vtu)‘^ + g^'"{Yru)‘^ + I+ V (r-)u^ ) r'^drdxdy, 


Fr{u-,[Ti,T2])-. = 






g'^ (VtU) + 9^^tu){Vru )) r^dtdxdy. 
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We now wish to take the limits Ri —?■ r_|_ and oo in the fluxes so that energy is 

defined across the whole exterior of the black hole. Approaching the event horizon we 
get 

lim Tr^{u][Ti,T 2 ]) = [ g''\Vtufr 1 dtdxdy=:F{u][Ti,T 2 ]), 

Mn,T2] 

we observe that is positive on the horizon making this flux a positive quantity. 

As for the contribution on X we can quickly see from the asymptotics that for Dirichlet 
and Neumann data that 

lim [ri,r 2 ]) = 0, 

H2^00 

however for time independent Robin data 


lim Tr^{u-, [ri,r 2 ]) = - lim 

R2—^oo r^oc 


r A 


'^u)^)dtdxdy 


1 

w 


[ (r2 ^u)'^l3dxdy 

Jt^ 

i 1,00 


1 


(r 


2 ^u) 


’(ddxdy, 


3 

where we understand the terms in the integral at r = 00 to mean lim.r_>.oo r^~^u. We 
may thus dehne the renormalised energy for a function u as 


Et[u] = \ j 2 rJ ^ r^drdxdy 

X,oo 


where we include the latter term for Robin data. This energy is positive definite and 
finite for our boundary data. It also satishes the useful identity 


(7) EtM=EtA^]-F{u-ATi,T2]), 

so that Et[u] is a strictly non-increasing function of t. Using the redshift effect [16], we 
can remove the degeneracy at the horizon to establish: 

Theorem 4.1. Suppose u is a solution to © satisfying suitable (Dirichlet, Neumann or 
Robin) boundary conditions at infinity. Define the non-degenerate renormalised energy 
density S[u\ by 

(8) £[u] := + r^(Vru)^ + iVtu)^ + r‘^\^u\^. 

Then there exists a constant C = C{M,1 ,k) > 0 such that for any Ti < T 2 we have: 



£\u]drdxdy < C 



£[u]drdxdy 


Remark 4.1. If enough regularity is assumed on the initial data of 0 one may extend 
this result: by applying T and the red shift vector field as commutators, together with 
elliptic estimates and a Sobolev embedding we can extract pointwise boundedness. In this 
setting one can prove a result similar to that found in |8| . 
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5. Decay Rates 

Now that we have established energy decay for solutions to Q we turn to the problem 
of establishing integrated decay. Due to limitations of our energy current we show the 
existence of a k* with | < k* < 1, such that our results hold for k G (0, k*). We believe 
the results hold for all k G (0,1) however the energy current will be far more complex. 
We will examine all the previously listed boundary conditions. We remark that this 
range of k includes the conformally coupled case (k = g). The methods we will use were 
first established by Morawetz [20] and m for the obstacle problem for waves and also 
seen in [12], [16] and [2T] for the Schwarzschild black hole. The core idea is to repeat a 
similar argument as in the energy decay proof, only rather than using dt as a multiplier 
we examine vector fields of the form h{r)dr for some radial function h. 


5.1. The Morawetz Estimate. In this section, we shall establish the following inte¬ 
grated decay estimate 


Theorem 5.1. [Precise Version] Suppose u is a smooth solution to ([^ with k G (0, k*), 
satisfying Dirichlet, Neumann or Robin boundary conditions at infinity. Then for any 
Ti < T 2 we have: 



dr] <C f 
J'Eti 


S[u\drdxdy 


with C = C{M, I, k) > 0. 


To prove this we make use of the following two lemmas 

Lemma 5.1 (Divergence of the Modified Energy Current). Let u be a solution to (§ 
with Dirichlet, Neumann or Robin boundary conditions. For an energy current defined 
as 

>[m] = + wi(r)uV^u + W 2 (r)u^X^, 

where 

X = rdr, 

we then have 


X^J^[u] = ( 2rw2{r) + 


r^w'fir) 2Mw'i{r) 


P r / r'^ 


uVrU + 


2pMw[{r) 


-uXtu 


+ ^W2{r) + 2KW2{r) + 


(3 — 2k)^M 2{3 — 2K)‘^Mwi{r)\ 2 




+ 


4^3 


u 


^ ((2 — k)PM PMwi{r)\ ^ ^ 1 ^ I 

(9) +A[ + VruVtu + ((1 - k) + wi{r)) \ fu\ 

\ O j 

Kr‘^ (3 — 2k)M r‘^wi{r) 2Mwi{r)' 


+ - 




{VrUY 


P r P r 

{5-2k)PM 2pMwi(r) \ 

+ 1 - —-r— - -y ((1 - k) + u;i(r)) {Vtuf. 


The proof of this may be found in the technical lemmas section of the paper. 
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Lemma 5.2 (Bounded Integrated Divergence). For a current defined as 

( 10 ) P[u] = + wi{r)uV^^u + W2{r)u‘^X^^, 

where 

• u solves with Dirichlet, Neumann or Robin data. 

• X = rdr, 

• Wi = -ki + f{r), 
with ki > 0 and f G 0{r~^), 

• W2 = ^, 

with 0 < k 2 < M. 

We have that 

f -Vf,P[u]d'^Vol < C 

for a constant C independent of Ti and T 2 . 

The proof of this may be found in the technical 


£[u]drdxdy 


-Ti 


lemmas section of the paper. 


Proof of Theorem 5.1 , We split the estimate into two parts as inspection of ([^ reveals 
that the coefficient of the tangential derivatives appears with the opposite sign to the 
higher order terms of the time derivative’s coefficient. This is problematic when trying 
to get a signed divergence. So we initially cancel this term off to control a positive 
quantity. This establishes part of the Morawetz estimate and is then used to control a 
divergence that includes poorly signed terms of lower order but crucially correctly signed 
tangential terms. 

We define our first current as 


with 


Jf [u] = + {n- l)uX^^u + 


X = rdr 


{3-2k)M 

2j.3 


1 e 

2 2 


u 


2xu 


as our multiplier. It can be easily checked that this current satisfies Lemma [5 .2 1 providing 
e < 2(1 — k) and by using the divergence formula from Lemma |5.1| we get 




- 




r 


XrUXtU + 
4 


(3 — InfiMe o 

-i-^ 

4r 


— -(3 — 2k)(1 + e)MuVrU + (Vru) 


which we can re-write as: 


( 11 ) 




= 


3l^M 


(Vtuf - 




VrUXtU + -prC^rUY 
r 


+ ^^ e)MuVrU + Mr{Xruf. 
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Note that our choice of tci ensures that there is no term involving ^u. We now seek to 
bound a positive quantity from above by this divergence. First, note that 


AfM 


VrUVtU 


< ^^(Vtn)2 + 5Mr(V,u)2, 
dr'^ 


for some <5 > 0 to be later determined. We thus have 
,2 4/2m~ _ rV- ,2 

- ^{VtUy - VrUVtU + -p^{VrUy > - 




P 


3 - T i^tuY +{^-5Mr] {Vruy. 


p 


This deals with the first line of 0- Now, for the second line 
(3 - 2K)2Me 2 1 


4r 

M f 


-u^ — ^(3 — 2 k )(1 + e)MuVrU + Mr{VruY 


3-2kV^2 /5/^3-2Ky 2 (l + e)2^2 


u — 


2/3 


r^^VrUf + r^{VrU)^ 


where /3 > 0 is a constant from Young’s inequality that we shall later determine. Fac¬ 
toring the above we deduce 

(3 - 2K)2Me 2 1 


4r 


-u^ — ^(3 — 2 k )(1 -|- e)MuVrU + Mr{VruY 


> 


M 


e — 


/3 


3 — 2k 
2 


+ 1 - 


(1 + 6)- 
2/3 


r‘^{Vru) 


now combining everything we have 


-V^JfM-r2> 

so providing that 
( 12 ) 
and 
(13) 


Pm 


M 


3--]{Vtuf+[^-Mr[5 + 


(1 + 6)- 
2/3 


- 1 


(VrU) 




3 — 2k 


/3<2e, 


u 






we have a bounded positive quantity. 

We now show for which range of e this can hold and then convert to the mass range this 
method works for. 

Firstly note (13) is equivalent to 


<5 < 3- 


(1 + 6)- 

2/3 


4 

3 


< 3 - 


(1 + 6)^ 

2/3 


but for consistency we then need 
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which is equivalent to 




again, for consistency we then also need 


which is equivalent to 
which holds for the range 


^(l + e)^<2e, 


3e^ - 14e + 3 < 0, 


G Q(7-2\/Io),^(7 + 2\/Io)) , 


now we also have the restraint that 


that is 


e < 2 - 2 k, 


-=<1-2, 


which is extremized for 




K^vro 


- 1 PS 0.887, 


thus providing an upper bound for the k ranges we can currently prove for. If we also 
make the choice 

e = 1 — K, 

we can easily then show the positivity and boundedness for the range 


K G ( 0, ^ (2^10 


- 1 


So for K in this range we may find all the constants to deduce that 


1 


-r^^c + r\Vruy + -u^j , 

for some c = c(M, I, k) > 0. We remark that we have used r > {2MP)3 . 
We now invoke Lemma 15.21 to deduce 


(14) 


IM 


ITi,T2] 


^{Vtuy + r^{Vruy + ) dtdrdxdy < C* / £[u]drdxdy. 

/st, 


Now that we have control of u, Vtu and in an integrated sense, we can return to 
establish an estimate involving ^u. To obtain this, we consider the current: 

J^ = -(\ + {l- kU uV>^u. 
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The divergence of this current is readily calculated to give: 


-V^4[u]-r^ = 
+ 


2 /6M 3\ ~ _ 4Z2M(r3-l) 

= I-I ii\/rU -P— UVtU -T- -VrUVtU 




P 


4 I 


P 


+ M r - ^ (V,u)2 + 


{Pm (3r3 - 2) - Pr^) 


{^tuf 


(3 - 2KfM P - 2) 
4r^ 


+ -\fu\ 


Examining the coefficients of the terms involving u, Vtu and we see that we already 
control all of these terms with appropriate weights. In particular, it is clear that we can 
find a C > 0 such that: 


- \fu\" < -r^ + C y-^{VtuY + r\Vruy + -u 

Integrating this estimate, applying Lemma |5.2| and estimate (0, we finally conclude 

+ r^{yruY + \{ytuY+ ^Wu\^\ drj <C I £[u]drdxdy 
\r r-i r J 


for some C > 0. 


□ 


5.2. Integrated Decay Estimate without weight loss. We can now restate our 
result regarding the integrated decay with no loss in the radial weights: 

Theorem 5.2. Suppose u is a smooth solution to Q with k G (0, n*), satisfying Dirich- 
let, Neumann or Robin boundary conditions at infinity. Let £[u] be the renormalised 
energy density of the field u, as in Then 


£[u]drj < C / {£\u] + £\ut\) drdxdy 


(15) / 

4>1[Ti.T2] 

for some C = C(M, I, k) > 0. 

In order to prove Theorem |5.2[ we shall require the following result: 

Lemma 5.3 (Hardy Estimate). Let 4> '■ [^+)C>o) — )• M 6e a smooth function such that 
limr^oof^(f) = 0. Then the following inequality holds 

'^dr<Cl f —dr + 


(16) 


{'Vr<f>)‘^r‘^dr ) . 


' r+ \J r+ ' JR 

Where C = C{R,r+) >0. 

The proof of this may be found in the technical lemmas section of the paper. 

Remark 5.1. For solutions to (§ we have from our earlier asymptotics for Robin data 
that r^u ^ O (r“^+'^) so the above Hardy estimate holds. 


Proof of Theorem \5.S\ In order to improve the weights in Theorem 5.1 we first apply 
the Hardy estimate of Lemma 5.3 to establish that 


(17) 


/ u'^dy < C £[u\drdxdy, 
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holds under the same conditions for u as in Theorem 5.1 Next we exploit some of the 
spacetime’s symmetry. We note that 


(18) 




= 0 , 


this tells us that (V^u) satisfies (provided the initial data is regular enough) and so 
0 applies. This yields an estimate of the form 


/ {'Vtu)‘^dr] S[ut]drdxdy. 

J AilTi,T2] 


we may then recombine with the estimate in Theorem |5.1| and see 


(19) / (Vtu)^ + r^(Vrtt)^ + u^dr] < C* / (T[m] + S[ut]) drdxdy 

it now only remains to recover the tangential derivatives. We will proceed with a more 
robust estimate here but refer the reader to the remark after this proof if they want a 
faster but less frugal route. 

We now define the current 

J^[u] = f^f^[u]X'^ - (2 - K)nV^ii 

where 

X = rdr 

from examining we can see that this current won’t pick up any cross terms. More 
explicitly we get 


. r" = 


l*M 


2r 


- ) {XtuY +{-^-Mr] {VruY + 


\2 I M{3 — 2kY 2 I 2 


4r 


-u 


+ r^\f 


u\ 


this current clearly satisfies the conditions of lemma |5.2| and as we control all the non- 
tangential terms we have that 

/ Y\^u\^dr]<C / {£[u] + £[ut]) drdxdy 


'^\TiX2\ 

which we combine to get 


•^A^[Ti,T2l 

and hnally 
( 20 ) 


iytuY +iyruY +u\^drj < C / {£[u] + £[ut]) drdxdy 


-Ti 


/ £[u]dr]<C {£[u] + £[ut]) drdxdy. 

4yVl[Tj,T2l 


□ 


Remark 5.2. If we aren’t concerned about the robustness of this result we can also 
exploit 

a' 

P\ 


( 21 ) 

( 22 ) 


dx, Dg + ,9 




a 

P\ 


= 0 , 
= 0 . 
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These would then manifest as VxU and VyU control. This would give us an estimate of 
the form 




{£[u] + £[ut] + £[ux] + £[uy\) drdxdy 


and this would be enough for a decay estimate. With nonlinear applications in mind, 
avoiding explicit use of the toric symmetry is preferable. 


In order to establish a decay statement, we require the following straightforward 
corollary of Gronwall’s Lemma; 

Lemma 5.4. Let k > 0, and suppose that f G (^^([rijoo)) satisfies 

(23) fit) < ->if{t) +- - -r, 

then there exists C = Cik, x) > 0 such that 

fit) < +- — -r- 

^ ^ ^ ^ il + t-Ti)^ 

The proof of this may be found in the technical lemmas section of the paper. 


We also require the following quantitative version of the redshift effect, which can be 
found in Theorem 3.8 of [22] . 


Lemma 5.5. There exists a modified energy the redshift energy such that 

• Et[u] is equivalent to the non-degenerate energy at time t. That is, for any smooth 
u satisfying appropriate boundary conditions we have: 

C~^ I £[u]drdxdy < Et[u] < C / £[u]drdxdy 

J-Bt 

for some C > 0. 

• If u solves S subject to Dirichlet, Neumann or Robin boundary conditions, then 
Et[u] satisfies 

(24) jEt[u]<-xEt[u] + CEt[u] 

for some x > 0. 


Theorem 5.3. Suppose u is a solution to ([^ with k G (0, k*), with Dirichlet, Neumann 
or Robin data. Then 


(25) 


f c 

£[u]drdxdy < 


/ £[d^u\drdxdy 


k=o • 


for some C = C(n, M, I, k) > 0. 


Proof. We can easily see that 


Et[u\ < C' / £[u]drdxdy 
JY.t 
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as the quantities are equivalent norms for u away from the horizon. Integrating in time 
we get 

f Es[u]ds < C / £[u]dr] 

JTi J AilTi,T2] 


SO by theorem 5.2 we have that 

f Es[u]ds < C [ {£[u\ + £[ut]) drdxdy. 

J J 

Since Et[u] is non-increasing, a slight adaptation of |23j . Lemma 5.8 immediately implies 
that for any t > Ti we have: 

Et[u] < ^ ^ rr [ + £[ut]) drdxdy. 

J'Zti 


I+ t-Ti 


This gives us decay of the degenerate energy. To establish decay of the full non¬ 
degenerate energy we apply Lemma 5.4 to (24), and we deduce that 

{£[u] + £{ut\) drdxdy. 




C 


l + t-Ti 


-Ti 


< 


-Ti 


£[u]drdxdy + — 


C 


+ t-Ti 


{£[u] + £[ut]) drdxdy 


-Ti 


i + t — Ti 
Finally, we obtain; 


{£[v\ -|- £[ut]) drdxdy 


C 


£[u]drdxdy < ^ ^ 


{£[u] + £[ut]) drdxdy. 


-Ti 


Which gives the result for re = 1 on setting Ti = 0. For higher re, the induction argument 
follows precisely as in [23], Lemma 5.8. □ 


5.3. Gaussian Beam and Derivative Loss. We now have an integrated decay esti¬ 
mate but with derivative loss. While we have not quantified precisely how much deriv¬ 
ative loss is required we will show that it is necessary. That is 

Theorem 5.4. There exists no constant C > Q, independent ofT, such that the estimate 
(26) j £[u]dri < C I £[u]drdxdy 

holds for all smooth solutions u of (§. 

The proof of this comes from the Gaussian beam construction as seen in [T0| . The core 
idea is to show that there exists null geodesics that can remain outside the event horizon 
for arbitrary lengths of coordinate time. We then can construct approximate solutions 
to (j^ (Gaussian beams) supported in a tubular neighbourhood of these geodesics such 
that they lose arbitrary small amounts of energy along them and remain close to true 
solutions in the energy norm. As we can find solutions which lose arbitrary small amounts 
of energy for any fixed time interval we cannot have an estimate of the form in Theorem 

[531 
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Lemma 5.6. For any given T > 0, there exists a null geodesic 7 such that Im{'y) C 
-^[o,r]n{|r+ < r < R{T)} for some large R{T). Furthermore 'y is a smooth embedding. 

The proof of this may be found in the technical lemmas section of the paper. 


We now use a modified theorem from m- These modifications are due to the fact 
that we are not working on a globally hyperbolic manifold but due to our geodesic being 
a smooth embedding an analogous result holds. 


Theorem 5.5 (Gaussian Beam). Let {M,g) be a time oriented Lorentzian manifold 
with time function t, foliated by the level sets = {t = r}. Furthermore, let y be a 
smooth geodesic embedding that intersects Sq and N a timelike, future directed vector 
field. 

For any neighbourhood M of y, for any T > 0 with S-p n Im{y) / 0 there exists a 
Gaussian beam u\ of the form u\{x) = with the following properties 

where the constant C(T) depends on ajy, (j) and T, but not on X, 

Et[u\] 00 for A —)• 00 , 


and 

u\ is supported in J\f, 
provided we have on Af [o,t] G H Sq), 

<6-, g{N,N)<-c<0, -g{N,nj:^{t))<C, 




and 


\g{Vns^N,ei)\,\g{Ve,N,ej)\ < C, for 1 < i,j < 3, 

where c and C are positive constants and 61 , 62 , 63 } is an orthonormal frame. 

In our case the vector field N is dt and the time function is simply the t coordinate. 
With the null geodesic from lemma 5.6 all the conditions can be easily seen to hold 
provided we bound Af away from I. 

We will also need the following result about (p from jlO] (2.14) 

Lemma 5.7. 3 '((/>|-y) = 3‘(V^|..j,) = 0 andXs{VV(p\ry) is positive definite on a 3-dimensional 
subspace transversal to y. 


Combining the fact that ajy is independent of A with (5.7) we have that the norm 
of our Gaussian beam is independent of A. We collect this observation in a lemma, 


Lemma 5.8. Let ux be the function constructed in 5.5 then there exists a constant 
C{T) > 0 independent of A such that the following bound holds: 

(27) < C(T). 

Lemma 5.9. For all e > 0 there exists a solution v of for all boundary data types, 
initial data supported away from the horizon, with = 1 and a Gaussian beam u\ 

such that 

\Et[v] - Et[ux]\ < €, VO < t < T. 
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Proof. Firstly construct u\ from theorem 5.5 using the geodesic in lemma 5.6 ensuring 
that Af is bounded away from I and Ti, and then define 

u\ 


u\ ■= 




so we have (by the triangle inequality) 


a 


□gMA + pUX 


L\M[o,t]) 


as A —)• oo. 

Now set V to be the solution to 


1-1 Q; 

UgV + —V = 0, 


v|so = ''^a|eo> 
nEo^l^o = '^So«a|eo> 

with Dirichlet, Robin or Neumann boundary conditions. 

We now apply the same energy estimates as in section 4 to ux, (we remark that as 
ux doesn’t solve (j^ this introduces an inhomogeneity in the bulk integral) yielding the 
inequality 

Et[ux] < Eo[ux] + (ngfiA + ^^a) VtUA 




after a simple application of Cauchy-Schwartz we hnd that 


Et[ux] < C{T) Eo[ux] + T ■ sup {Et[ux]) 

V ie[o,T] 


a 


OgUX + pUx 




for all 0 < t < T. (Note that Af being bounded away from infinity controls the flux 
there). 

Taking supremums on both sides and absorbing an application of Young’s inequality we 
get 

Et[ux] < C{T) (eo[ux] + \JgUx + ^ux 
Applying this inequality to the difference v — ux, gives us the result 

Oi 

\Et[v-ux]\<C{T) Ugiix + rp^ux V0<t<r, 

as they agree on Sq. So for fixed e we simply choose Aq large enough and set ft := uxq 
and V = fAp. □ 

We now invoke Theorem 2.36 from m which tells us that the Gaussian beam energy 
is localised around geodesic energy. That is 

Lemma 5.10. For all e > 0 there exists a neighbourhood A/q of N such that 

\Et[uxWo] - {-9{T,j)\im('r)nEt) I < 


for all 0 <t <T. 
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Lemma 5.11. Let T > 0 then for all e > 0 there exists a solution u of ^ whose initial 
data is supported away from the horizon, that can satisfy any choice of the discussed 
boundary conditions and with Eq[u] = 1 such that 

Et[u] > 1 - e, 


for all 0 < t < T. 


Proof. As T is Killing we have that —5(7", is constant. We may choose this 

constant to be 1 when we solve for 7 (see the proof of lemma 5.6 in the technical lemmas 
section for details). Applying the triangle inequality to the results in lemmas 5.9 and 
5. 10| yields the result. □ 


Proof of Theorem 5.f. We proceed by contradiction and assume there exists a constant 
C independent of T and u such that 

rT 


(28) 


Et[u]dt < CEq[u], 


as the energy is decreasing we have 

f Exiujdt < f Et[u]dt < CEq[u 
Jo Jo 

this gives us the estimate 

TEt[u] < CEo[u], 

now let T = 2C and construct u from lemma 


5.11 


where we choose e = | we then have 


5C<C. 


which is clearly a contradiction. Now that we have established that (28) cannot hold we 
may extend the result to the non-degenerate energy. 

If u is supported away from the horizon then there exists C > 0 such that 

£[u] < CEo[u], 


'So 


and clearly 


Et[u]dt < C 


IM 


£[u]dr}. 


[0,T] 


Assume by contradiction that the statement of Theorem 5.4 holds true. Then the above 
inequalities entail equation (28) contradicting the hrst part of the proof. □ 


6. Technical Lemmas 

6.1. Proof of Lemma 15.31 

Proof. Firstly define a cut-off function 


X{r) 


0 if r < i?, 

1 if r > 2R, 

Smooth if r G [R, 2R]. 
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With the property x^(r) < — for some C > 0 and monotone. 

This can be achieved by a suitable bump function. 

Now we write 

\W\\l^ = \\xu + {I - x)u\\l2, 

After applying the triangle inequality we estimate the terms separately 


II( 1 -x)«IIl 2 = / {l-x)udr< 
Jr+ 

And for the other term 


'•2R y2 

r • —dr < 2R I —dr. 


roc ^2 


’r+ 


'r+ 


/ oo noo noc 3 \2 / j.-2+2k\ 

ixufdr = ixufdr = i^xur^-^j dr 1^ 2^- 2 ) 


1 

2k-2 


ixufr 


+ Z 


J i? ^ JR 


/■“ - 1 

/ Xu{VrXu)rdr< - - \\Xu\\L2(r+,oo)\\r^rXu\\L‘^{R,oo)- 

JR L — K 


=0 


Now looking at 


Ir^rXuWhf^R^oo) = / r'^i^rxufdr 


IR 


< C 


roc 

/ r^ixVruf 
JR 


+ r‘^{udrx)‘^dr 


< C 


IR 


3 _ V? 

CVru)‘^r'^dr + / —dr 

Jr+ r 


Combining all these estimates yields (16). 

6.2. Proof of Lemma 15.11 

Proof. 

We recall that 

u\ X 


+ X-'S,. 


flU ■* 


We start with the latter term. A quick computation shows us 


1 


2 I (3 


V nUV^U 


and thus 


S'r = - { (3 — K)V{r)u'^ + 


X'^S^ = (3 - K)V{r)u^ + 


□ 


We now deal with the deformation tensor term. We can factor the metric out of the 
deformation tensor 


X 


TT = g 


3M , 2 SPm , , P(5pM + r^) , 2 
- dr - T^dtdr -^-E-idrC 


Contracting g into T we get 


= -V^nV^u - 2V{r)u^, 
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then combining with the S terms we get 


= (1 - K)Viry + 


Contracting the remaining deformation terms we get 




P {sPM + r3) 


iPM + 


^ruvtu - - 


2 M(3 2k)^ 2 , l|y*„.|2 


2^5 


2/2r 


4r^ 


So for the first term we conclude 


=- 


(2k — 3)^M 2 _i f ~ 3)M Kr'^ 


^^3 


- X I 


MW, |2 4:{k-2)PM~ ^ 

+ (1 — K)\^Uf - 5 - VrUVtU + 




Now we compute the divergence of the second term. To do this we write 


V,,u = I VnU — uV„l 


11^ , 


then we can compute 

[wi(r 


= V fj,(wiu)'V^u + wiuV fjV^u, 


now we explore part of the second term 


VfjX'^u = OgU + Vfj, (^uV^l 

= UgU + uV^V^l + V^uV^l 

= \I\gU + uVg + VgUg^’^Vul 


= OgU + nV. 




— + X- « —VaU 


fir 




= \JgU + - - At uVn - + X - - V„U - uVnl 


9‘ 


r 

fir 


= OaU +- k] uV 






2 TV 

g 


u + 


(i-K) 


a 3 (I -At) 2M{l-Kf (i-At)^ 1~ 

= -^u + —’-u + -- ’—u - ’ u + -V^u 


V^u 


P ' P 


p 


r 
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Recombining we get 

f2M 1 ~ 

V^{wiu)V^u + wiuV = (Vrici) uV^u + wiV^uV^u + wiu ( —H—V^’u 

= {VrWl) uV^u + WlV i_iuV^u — wiuV i_ilV^U 


+ WlU 


/2M 

( 7 ^ 


u + - V'’m 


2M 


= < 7 ^'’ (VrWi) uVfj^u + wig^’^VnuVlyU + wi — 


Explicitly this is 

2pMw[{r] 


uVtu + 


r‘^w\{r) 2Mw\(r)\ ~ ApMwi(r)~ ^ 


P 


+ 


r‘^wi{r) 2Mwi{r) 
P r 


r J 

;2„ 


(V,n)^ + - 


2rMwi{r) rwi{r) 


(Vtu)' 


/ n2M 2 , MW, ,2 

+ wi{r)^^u + w\{r)\fu\ . 


For the last term 

V^{w2{r)u^X^^) = \I^{W2U^)X>^ + W2V?X^X^^ 

= X f_i{w2)u^X^ + 2w2UV^uX^ + 3W2u‘^ 

= {rw2{r) + 3w2{r)) + 2rw2{r)vX/rU 

= {rw' 2 {r) + 2 KW 2 {r)) v? + 2 rw 2 {r)uXrU 

Combining all of these results and tidying up the algebra we may now fully express the 
divergence for the current as 

- ( r‘^w\ir) 2Mw\{r)\ ~ 2pMw\(r) 

X^J^[u] = \ 2rw2{r) H-- - -j H-- -uXfU 


+ ^W2{r) + 2KW2{r) + 


{3-2KfM 2{3-2KfMwi{r)\ 2 




+ 


^^3 


(29) 


+ 4 


+ - 


{2 — k)PM PMwi{r)\ ~ ^ iw, i 2 

- 7 -^- 7 - j ^ruVtu + ((1 - k ) + Wi{r)) 

{3 — 2k)M r‘^wi{r) 2Mw\{r) 


P 


+ 


P 


{Xruf 


+ - 


{5-2k)PM 2pMwi{r) I 


- Z9 ((l-^)+^i(^)) i^tu) . 


□ 


6.3. Proof of Lemma 15.21 

Proof. To prove this we apply the divergence theorem (|^ to the current and prove the 
fluxes across the hypersurfaces are bounded by a constant multiple of the initial energy. 
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For the surfaces of constant r we compute that: 


(30) 




(3 - 2k)2 
4 



^ ^ 2P 


2Mf) {Vruf + 




— + wi(r)^ (r^ — 2Ml^) uVrU + uVtu dtdxdy. 

2 ^ r 


(Vtu)2 


We first investigate the terms on X. Using the asymptotic analysis for Robin data as 
seen in Q we can quickly see most of these terms are converging to 0. The limit thus 
reduces to the study of 



2M/2) (Vru)^ + wi{r)^ (r^ - 2Mf) uVrU 



where in the last line we have substituted in the limit for Robin data and A > 0 is some 
constant independent of T. In integral form this is 


lim / 

r^oa Jy;IXi,T2I 


J^m^dSY,r 


lim / 

1—>■00 


xp 



2 


X j f /3 (r^ dxdydt. 

JTl 


As for the contribution of the terms at the horizon (r 
Young’s inequality 


r+) we first bound (30) using 


lim / 


lim / 


Jum^dSj:^ < 


(3 - 2k) 


I \ 2 

-M + e j ^ (r^ - 2Mf ) {Vruf 


r'^W2{r) - ^ 

+ ( —5- + ^ H-5-) (Vtu)^ - -r^\^u\^ + wi{r)-z (r^ - 2 MP) uVrudtdxdy 

\ 2 er^ J 2 ' 


< 


qTi.Tal 


r^W2{r+) - 


(3 - 2k)^ 


M + e W + 


Pm Pr+ PM‘^wl{r+)\ 2 


2 “I-A-^ 

rj' 2 


er^ J 


tuy dtdxdy. 


< 




k 2 - + e]u^ + CFu[Ti,T 2 ] 


For k 2 < jVf -^e can always find an e > 0 such that the first bracket in the 

integrand is negative. We may drop the negative terms and evaluate all the radial 
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functions absorbing their values into a constant C. We thus conclude 






C{Vtufdtdxdy < CFu[Ti,T 2 ] < CEtM- 


We now have expressions for the X and contributions. 

As the energy is decreasing in time both surfaces of constant t can be dealt with in one 
calculation. So for Ti < t < T 2 we compute 




= / -2fM{Vruf + 


2l^M 


+ l^r VruVtu + wi{r) 


2l^M 


+ r ) uVtu 


2pMwi{r] 


uVrU drdxdy 


2l^M \ ~ 

= J —2pM(Vru)^ + —2 - I'^r j VruVtu — kil'^uVtu 

f2fM{f{r)-ki) , 2 ,, ^ 2pMwi(r) ~ , , , 

+ ( - 5 - - + f‘f{r) ) vSJtu - -uVrU drdxdy 

\ r'^ J r 

= / —21^Miyru)^ + hi{r)VruVtu — kil^vSJtu + h 2 {r)vSJtu — h^{r)uVrU drdxdy, 
J-Et 

where hi{r) G 0{r), h 2 {r) G 0(r~^) and /i 3 (r) G 0{r~^). With this in mind we now 
apply Young’s inequality to the cross terms 

[ J^m^^dSE.t 


+ ^ + Gl^ 2 (r)|^ {Vtuf + (l\h 2 {r)\ + 7 r'^|/r 3 (r’)ll drdxdy 


2ei 2e3 2 


kfl'^es 


+ / y-j\h 3 {r)\ +—{hi{r)) -21 Mj (Vru) H- - — u drdxdy. 

J Sj 

Now r is bounded away from zero and as the function coefficients are of a low enough 
order (in terms of r) we can easily bound the first integral by a constant multiple of the 
energy at time Ti (the constant depends on the choice of the Cj’s). To deal however with 
the second integral we need to invoke the Hardy inequality 

[ Jam^dSE,, 


-Je ^ ^ ^ ^l^2(r')|^ (Vtuf + Q|/i2(r)| + ^l^ 3 (r)|^ drdxdy 

+ J + ^\h 3 {r)\ + ^(/ii(r))^ - 2l^M^ {Vruf drdxdy. 

We now note that we simply need to choose positive ei, €2 and 63 such that 

(j2h^rl + ^|/i3(r+)| + |^(hi(r+))2 < 2fM, 
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rT2 


we may then find a C > 0 independent of Ti and T 2 such that 

[ J^m^dSY.^ < CEt^[u], 

J-Et 

for Ti < t < T 2 . 

We now apply the divergence theorem to yield 
[ -V^J^^[u]d‘^Vol < CiEt, [u]+C 2 Et, [w] +C 3 ET, [u ]-C 

which may be rewritten as 

f j I3{ru)‘^dxdydt + f — < C-Eti M, 

thus proving the proposition. 




/3{ru)‘^dxdydt, 


□ 


6.4. Proof of Lemma 15.41 


Proof. We first re-write (23) as: 

fL 

dt 


so that 


Here, we use: 


d Ap^ 

I [f{t)e^) < 


{l + t-TiY 


Jt 


-ds 


Ti (l + S-Ti)^ 
rt-Ti Ks' 

/ -rds' 


< CAe.^'^^ 


0 {l + s'Y 

/ gx(t-Ti) 


A 


lo (1 + s')' 


-ds' = 


>c{i + s'y 


ft pXs' 

+ k / — -fords' 

Jo X (1 -I- s') 




< - p + kt max 

x(l-bt) ^'6(04) 

pXt 

< c- 


X (1 -h s') 


fc+l 


( i+tr 


a simple re-arrangement gives 


□ 

















6.5. Proof of Lemma 15.61 


Proof. The easiest way to construct the geodesic is by the Hamiltonian method. We can 
quickly spot three integrals of motion arising from the Killing fields dt,dx, dy and from 
the fact that 7 is null. More explicitly we define three constants a, b, c and write the 
equations 

a = g{i,dt), 
b = g{i,dx), 
c = gii,dy), 

0 = 5 ( 7 , 7 )- 


Working in co-ordinates and taking an affine parameter r for 'y(t,r,x,y). We get the 
following geodesic equation 



where df = b‘^ + (? and we chose the negative root of as we want to look at a photon 
travelling tangentially at a distance R from the origin falling into the black hole. 

r( 0 ) = R, 
r( 0 ) = 0 . 


This allows us to directly solve for a as 


a = —d 



2 

it) 


We have chosen the negative root as we would like a positive killing energy and i > 0 
so we can measure the co-ordinate time it takes to fall some distance towards the event 
horizon. Our equation thus becomes 


t = 



2 M/2 / I 1 \ 2 


( 32 ) 


r = —dV2M 
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So long as i? > r_|_ we have that these are all signed (away from R) quantities. This 
allows to deduce that t, r, x, y are all monotonic. As r is monotonic decreasing the geo¬ 
desic never reaches I and we find that t, x, y also are monotonic. As 7 is clearly smooth 
we deduce that 7 is a smooth embedding. 


We now wish to study the co-ordinate time it takes to fall a distance of ^ remaining 
entirely outside the event horizon. More explicitly let |r_|_ < ^ we are interested in 

Ant = t\r=R — t\ R, 

2 2 

To evaluate this we first need the affine time tr between R and =?. From (|32|) we can 

2 ^ — 

compute this as 


^2 d\/2M 


rR 


r3 ji3 ) 


rescaling the integral we get 


Now 


0 < 



Y- < 00 , 


so this integral is simply some positive constant. We denote it as Kp and write 


KpR^ 

Tr = - - - 

2 dy/AM 


We now turn to the fall time. Our goal is to show we can for fixed T, find a geodesic 
that remains outside the event horizon up to and including T. In order to do this we 
consider bounding i below on the interval r G ["f j -R] 


. dP 

t > -w 


2M 

~w 


> 


dP 


R^ - 2MP \ 5 

W / 


MRP I WMP 

{R^ - 16MP) ^ 

AdP 


R2{R3 - 16M^2)' 


t > 


now fix e > 1 and set C = I ^(l — y) then provided R^ > 2eMP we have 

dPC AdP 

rI{R^-16MP) 

C ^ V 

_ 16 M/ 2 ) J 

If we then choose e = 9 so we require R^ > 18MP we have 

dPl'^_ A \ 

V9^ RhMP 



r > 


)’ 
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Insisting > max{18MZ ^ lewl we find 

. Ml 1 

^ - 'mW' 


We may now prove an inequality for the fall time (for large R) 


Now fix T and let R 


f-TR . (IP 1 

Aflt = / ^ idr > K-^ ■ tr = KpR^ ■ 

2 Jq R 2 

max{|^, (18M/^)3, then have 


ArI > 

2 


3T 


or in other words a geodesic with r > |r+ for 0 < t < T. 


□ 
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